Quantum control of optical four-wave mixing with femtosecond u-3uj laser pulses: 

coherent ac Stark nonlinear spectroscopy 
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The four-wave mixing produced with two ultrashort phase-locked lu-3u> laser pulses propagating 
coherently in a two-level system in the infrared spectral region is shown to depend on the pulses 
relative phase. The Maxwell-Bloch equations are solved beyond the rotating-wave approximation to 
account for field frequencies which are largely detuned from the atomic resonance. The relative phase 
dominating the efficiency of the coupling to the 5lo anti-Stokes Raman component is determined 
by sign of the total ac Stark shift induced in the system, in such a way that the phase influence 
disappears precisely where the ac Stark effect due to both pulses is compensated. This fundamental 
quantum interference effect can be the basis for nonlinear ultrafast optical spectroscopy techniques. 
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Phenomena arising from the coherent control of non- 
linear interactions are of importance in fields as diverse as 
optoelectronics and materials research, in high harmonic 
generation, in photoionization or molecular dissociation, 
and in biological applications such as spectroscopy and 
imaging, among others 0, 0, El 0, EL • Optical quan- 
tum coherent control is based on the fact that the phases 
of interfering transition amplitudes in light-matter inter- 
actions can be controlled through the optical phase of 
coherent light sources that drive the interaction, in such 
a way that the transition rates to final states and the dy- 
namics at various stages of the process can be modified 



In a recent paper gj, a theoretical investigation on 
the quantum coherent control of the optical transient 
four-wave mixing of two intense phase-locked femtosec- 
ond laser pulses of central angular frequencies uj and 3lu 
propagating in a two-level atom (TLA) was reported. It 
was shown how the nonlinear (x^ ) coupling to the anti- 
Stokes Raman field at frequency 5lu depends critically on 
the initial relative phase <f> of the propagating pulses. In 
Ref. @, the study was centered to intense pulses in the 
visible and ultraviolet spectral regions, with frequencies 
at resonance or lower than the atomic transition. The 
phenomena observed in however, can be scaled to var- 
ious laser and material parameters. In the infrared spec- 
tral range, for instance, experiments on ultrafast molec- 
ular dynamics arc frequently performed by help of two- 
color pump probe nonlinear spectroscopy techniques. In 
this type of studies, due to the high intensities inherent to 
ultrashort (subpicosecond) pulses, nonlinear effects such 
as stimulated Raman processes may become important 
and are often utilized as a complementary tool to gain 
information [{|. 

In this Letter, we address the conditions for quan- 
tum coherent control of transient four-wave mixing in- 
teractions with two phase-locked uo-iu femtosecond laser 
pulses in the mid-IR spectral region. Our purpose is to 
reveal the basic physics for a strict two-level medium and 
to this end the Maxwell-Bloch TLA will be considered. 
We will study pulses with a duration of 300 fs (with a 
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FIG. 1: (Color online) Schematic energy level diagram. The 
resonance wavelength is considered in the mid-IR at A(^ 12 ) = 
3000 nm 



spectral width of 35 cm -1 ) and peak intensities as 
~10 8 W/cm 2 , which are typically used in infrared non- 
linear spectroscopy experiments |9( . We will examine the 
influence of frequency detuning by considering the field at 
3lu above resonance with respect to the atomic transition, 
something that was not considered in ||. Under these 
conditions, we will show that the quantum interferences 
leading to a dominating efficiency of the nonlinear cou- 
pling to the anti-Stokes Raman component are governed 
by the ac Stark shift induced in the system. Furthermore, 
we will demonstrate that the relative spectral amplitude 
of the anti-Stokes fields produced by phase-locked pulses 
cancels for frequency detunings that compensate the ac 
Stark effect. This fundamental interference effect might 
be the basis for a novel ultrafast spectroscopy tool based 
on coherent control, which we name coherent ac Stark 
nonlinear spectroscopy (CSNS) for use later below. 

The pulse propagation is modeled by means of the 
Maxwell-Bloch equations beyond the rotating-wave ap- 
proximation, allowing the resonant as well as the non- 
resonant regimes of the system to be described [l| . The 
equations are written as 



dH 
~dt 



1 dE 

Mo dz ' 



2 



dE 
~dt 
dpi 
dt 
dp 2 

at 
at 



1 dH N at n 
e dz e Q T 2 
1 

% 
1 



J2 



(1) 



r 2 
1 

"rl 



P2 



2 1 1 



Ep 3 - W12P1, 



(P3 - P30) 



2/y 



£,°2, 



where H(z,t) and E(z,t) represent the magnetic and 
electric fields propagating along the z direction, respec- 
tively, po an d eo are the magnetic permeability and elec- 
tric permittivity of free space, respectively, N a t = 2x 10 24 
m -3 is the density of polarizable atoms, p, — 4.2 x 10~ 29 
Cm is the effective dipole coupling coefficient, T\ = T 2 = 
1 ps are the excited-state lifetime and dephasing time, 
respectively, p\ and p 2 are the real and imaginary com- 
ponents of the polarization, and iO\ 2 is the transition res- 
onance angular frequency of the two level medium, con- 
sidered in the present simulations in the mid-IR region 
at 3000 nm (see Fig. P). The population difference is pa, 
and P30 represents its initial value. An hyperbolic secant 
two-color pulse that can be expressed as 



E(t)=E aJ (t)+E 3u) (t)=E sech((t-t )/t p )x 
[cos(uj(t — to)) + cos{Zuj{t — to) + ([>)] 



(2) 



is externally injected to the system. The peak input elec- 
tric field amplitude Eq is chosen the same for both pulses 
and results in an intensity of 4.0 x 10 s W/cm 2 . The dura- 
tion of the pulses is given by t p = r p /1.763, with t p = 300 
fs being the full width at half maximum (FWHM) of the 
pulse intensity envelope. t gives the offset position of 
the pulse center at t = 0, and it is the reference value 
for the phase of the pulses. The central angular frequen- 
cies of the pulses are ui and 3w, and cf> is the relative 
phase. The propagating system has been resolved nu- 
merically by means of a standard finite difference time 
domain method described elsewhere || . 

Figure 2 shows the field spectra at different propaga- 
tion lengths in the case that the central angular frequency 
of the field E^ is at resonance with the atomic transition 
(3u> = W12). The spectrum on the top is for the initial 
pulses. The succeeding plots show the evolution of the 
spectrum as the pulses propagate through the medium. 
In this case, we can observe the effect of the absorption 
of the pulse at 3a; , together with the appearance of other 
spectral contributions at 5u>, 7oj, and 9u>, which are pro- 
duced as a result of the coupling of the fields through the 
third order nonlincarity (x^) of the medium. It is clear 
that the conversion to the anti-Stokes Raman component 
(5u>) depends on the relative phase between the pulses. 
For the parameter values of the results shown in Fig. 
the coupling to the fifth harmonic (5w) is more efficient 
for (f> — than for <p — w, an scenario that was already 
reported in Ref. Q for a transition in the visible region. 



3. 

CO 

o 

CD 
Q_ 
Crt 

CD 




z=0 



z=10 nm 



z=25 nm 



z=50 nm 



A 1 



z=100 nm 



5co 



7co 9co 



FIG. 2: (Color online) Spectra of the total field at different 
propagation lengths as indicated. In the case shown, the cen- 
tral pulse frequency 3w is in resonance with the atomic tran- 
sition. Each spectrum is plotted in logarithmic scale. The 
relative spectral amplitude of the anti-Stokes 5u> components 
is shown quantitatively for z — 25 /im in Fig. |S| (case with 
77 = 9.0). 



Hence here we confirm that observed in Ref. |g for a dif- 
ferent system, particularly, considering a mid-IR atomic 
transition. It is important to note that this phase depen- 
dence effect involves the anti-Stokes Raman component 
(5cj) only, not the 7lo nor the 9w spectral components, 
which remain insensitive to the initial relative phase of 
the pulses. 

We now turn to the study of the frequency detuning 
of the fields with respect to the atomic transition. Wc 
will show that there is a central pulse frequency lu at 
which the relative phase dependence of the coupling to 
the anti-Stokes Raman component, which has been pre- 
viously discussed in Fig. [21 disappears. We will observe 
that this effect can occur because the ac Stark shifts pro- 
duced in the medium by the fields E u and E^ can be 
compensated when u> < lo\ 2 < 3a;. We will then conclude 
that the ac Stark shift in the medium governs the relative 
phase dependence of the coherent four-wave coupling. 

Indeed, the ac Stark frequency shift Ao; produced by 
a field of frequency to which is not near resonance in a 
transition of frequency u>i 2 can be expressed as [TT| : 
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where + is for u> < lj\ 2 and — is for u> > Wi2) (3 = 
p\E\/(2h), with \E\ being the amplitude of the field, 
a;i2 the frequency of the atomic transition, and lu the 
field angular frequency. The conditions f3/u> << 1 and 
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FIG. 3: (Color online) Spectra for differents values of r\ showing the anti-Stokes 5ui frequency component of the field at the 
propagation distance z — 25 ^m, for = (solid lines) and <fi = ir (dashed lines). At r\ = 5 (center plot) the dependence on 
the relative phase disappears. 



/3/o;i2 << 1 must be met for Eq. El to be valid Far 
from resonance, where 
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Eq. (PI becomes 



Ao;«/3 2 



10\2 - W 



^12 



(4) 



(5) 



Note that the last term inside the brackets in Eq. (JSJ is 
important far from resonance, where the rotating wave 
approximation does not apply. Clearly from Eq. 
when 10 > ujyi we have Aw < 0, and the separation in 
frequency of the states (wi2 + Aw) appears to be less than 
in the absence of the field. Contrarily, for to < lo\2 , the 
ac Stark frequency shift is positive. 

In the case of the present investigation, the main con- 
tributions to the ac Stark effect come from the two com- 
ponents of the propagating pulses E u and E 3u> , which 
have central angular frequencies to and 3a;, respectively, 
and the conditions for Eq. |JS} are clearly met for the 
parameter values of our study. Therefore, requiring that 
the combined ac Stark effect is null 
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we obtain to = co\2l 'v5- 

We will next show that the ac Stark cancellation fre- 
quency to — 1012/V5 is indeed observed with accuracy 
from numerical simulations. In Figure the spectral 
amplitude of the anti-Stokes 5a; component is shown 
for different values of the detuning between the fields 
and the atomic transition at a propagation length such 
as z = 25/im. It is useful to define the parameter 
7] = (uj^/lj) 2 , which sets the value of the detuning of the 
E u and E 3l0 fields. The spectra corresponding to <p = 
is shown by solid lines in Fig. while the spectra for 
4> = 7r is represented by the dashed lines. Clearly, there 
is a switch in the tendency to dominate the conversion to 



the 5a; anti-Stokes component. For 7/ < 5, the ac Stark 
shift due to E u dominates over the shift induced by E 3u] , 
and therefore the resulting separation of the states due 
to the combined Stark effect appears to be larger than in 
the absence of fields. In this situation, the coupling to the 
buj anti-Stokes field is enhanced for 4> = n, as it can be 
observed in Fig. |3 left plots. Contrarily, for 77 > 5, when 
the ac Stark shift due to E 3uJ dominates over the shift 
induced by E u , the coupling to the 5a; anti-Stokes Ra- 
man component is enhanced for = (right plots in Fig. 
|3J). In this last case the resulting separation in frequency 
of the states due to the total Stark effect is less than in 
the absence of fields. Furthermore, as shown in Fig. |3 
the switching occurs at r\ = 5, where the total ac Stark 
shift is cancelled, as expected from Eq. |SJ Note that for 
the numerical simulations to agree with the prediction 
of the analytical theory, the expression for the Stark ef- 
fect beyond the rotating wave ap prox imation needs to be 
considered (see Eq. (33) in Ref. [ll|h 

We now look at the relative spectral amplitude of the 
anti-stokes Raman fields, which we define as 



(7) 



with E^ UJ (z) being the spectral amplitude at the gen- 
erated anti-Stokes (5a;) frequency in the case that the 
initial relative phase is <j) = (see Fig. |3J), and E^{z) 
being the spectral amplitude at 5a; in the case that the 
initial relative phase between the pulses is <j) — n. Figure 
21 shows the results obtained from the numerical simula- 
tions for z — 25/im. Relative differences in amplitudes as 
|A_E5 W | > 0.5 can readily be produced in the cases con- 
sidered in our simulations. Moreover, we have checked 
that the relative phase dependence at rj = 5.0 remains 
null for propagation distances as long as z = 100/im. 

We have therefore demonstrated that considering the 
coherent propagation of two-color phase-locked femtosec- 
ond pulses in a two-level medium, with central angu- 
lar frequencies to and 3a;, one can find an angular fre- 
quency io at which the ac Stark effect produced by the 
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FIG. 4: (Color online) Relative spectral amplitude of the 
anti-stokes Raman fields [as defined in Eq. iJTJ] as a function 
of the parameter r\ for z = 25/^m. The vertical dashed line is 
■q — 5 as indicated. The dashed curve is a guide to the eye. 



propagating pulses is cancelled. At this frequency value, 
which requires the pulse with central frequency 3lu to be 
above resonance with respect to the atomic transition, 
the phase dependence of the transient four- wave coupling 
through the nonlinear susceptibility of the medium 
disappears. This is a fundamental quantum interference 
effect that to the best of our knowledge has not been 
reported before. 

In the present Letter, we have considered atomic fre- 
quencies which lie in the mid-IR region, which are of 
interest for applications e.g. in infrared nonlinear spec- 
troscopic techniques. It has to be stressed however that 
the phenomena that we report can be scaled to several 
material and pulse parameters. We can hence imagine 
straightforward applications for nonlinear ultrafast spec- 
troscopy techniques based on the coherent control of sub- 
picosecond two-color w-3o> propagating laser pulses. In- 
deed, the production of phase related w-3w pulses is fre- 
quently accomplished by some frequency trippling mech- 
anism with the subsequent variation of the phase of one 
of the pulses in order to obtain experimental control over 
the relative phase. Although our analysis has obviously 
been simplified by considering two well isolated levels as 
a first approach, the switching effect discussed here is of 
a fundamental level, and in that sense it should be ob- 
served experimentally in particular media where the two- 
level approximation is met. For instance, some gaseous 
atoms have well isolated resonances (as e.g. rubidium 
in the visible region). Also, the two-level approximation 
can be used for studying coherent effects in materials 



with a broad distribution of transitions, such as inho- 
mogeneously broadened resonance lines in gases and in 
condensed matter, and even for inhomogeneousquasi- 
continuous energy bands as in semiconductors |l2L Il8j | . 
CSNS can hence provide information on the transitions 
being probed by measuring the relative spectral ampli- 
tude of the anti-Stokes Raman fields with a simple scan 
of the laser frequency. The ac Stark mediated coherent 
control scenario reported in this Letter has therefore a 
broad interest and may be the basis for future studies on 
more complex systems. 
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